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Due to underlying symmetries the ground states of magnetic adatoms may be highly stable, which 
opens perspectives for application as single-atom memory. A specific example is a single holmium 
atom (with J = 8) on a platinum (111) surface for which exceptionally long lifetimes were observed 
in recent scanning tunneling microscopy studies. For control and read-out the atom must be coupled 
to electronic contacts. Hence the spin dynamics of the system is governed by a quantum master 
equation. Our analysis shows that in general it cannot be reduced to a classical master equation 
in the basis of the unperturbed crystal-field Hamiltonian. Rather, depending on parameters and 
control fields, “environment induced superselection” principles choose the appropriate set of basis 
states, which in turn determines the specific relaxation channels and lifetimes. Our simulations 
suggest that in ideal situations the lifetimes should be even longer than observed in the experiment. 

We, therefore, investigate the influence of various perturbations. We also study the initialization 
process of the state of the Ho atom by applied voltage pulses and conclude that fast, high fidelity 
preparation, on a 100 ns timescale, should be possible. 

PACS numbers: 72.25.-b, 05.30.-d, 03.65.Yz 


I. INTRODUCTION 

An ultimate boundary of miniaturization of informa¬ 
tion technology is reached when single atoms are used as 
memory bits. In this respect the experiment of Miya- 
machi et al. [l| with a single Ho atom positioned on 
a Pt(lll) surface (as depicted in Fig. [T|) represent an 
important milestone, promising lifetimes of several min¬ 
utes for the two degenerate magnetic ground states. For 
the crystal-field Hamiltonian with parameters as deter¬ 
mined in Ref. [I|, the two ground states of the Ho atom 
have ( Jz) ~ ±8 pointing into or out of the metal surface. 
The long lifetimes result from a combination of symme¬ 
tries and specihc properties of the system [l|. They are 
many orders of magnitude longer than those measured 
for single Co/Fe atoms on Pt(lll) or other mag¬ 

netic adatoms on a metallic surface [^. A related giant 
magnetic anisotropy of single adatoms had been observed 
earlier for Co atoms on Pt(lll) We should mention 
that in another recent work with Ho on Pt(lll), Do- 
nati et al. arrived at different crystal-field parame¬ 
ters, leading to ground states with (A) « ±6 and much 
shorter lifetimes. 

In this paper we develop the theory how to treat mag¬ 
netic adatoms on metallic surfaces which - for read-out 
and control - are coupled to electronic reservoirs formed 
by the metal substrate and a spin-polarized STM tip. 
The description of this dissipative system is based on a 
quantum master equation for the reduced density matrix 
for the angular momentum states, which accounts for the 
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FIG. 1. Setting of the experiment in Ref. [I|]: The magnetic 
state of an holmium atom on a Pt(lll) surface is studied by 
a scanning tunneling microscope with a spin-polarized tip. 


influence of the electronic contacts and the applied volt¬ 
age. For low voltages the description reduces to a rate 
equation in the eigenbasis of the unperturbed crystal- 
field Hamiltonian, which are specific superpositions of 
angular momentum states. However, for larger voltages 
the coherence of these superposition states is destroyed, 
and the relevant states are the steady-state eigenstates of 
the full dissipative quantum master equation. This is an 
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example for what has been called “environment-induced 
superselection” (or short “einselection”) principle [3, [1| ■ 
In the considered system the choice of the appropriate 
basis depends strongly on the voltage. It strongly influ¬ 
ences the relaxation channels and times. Switching the 
voltages also allows initializing the systems in a specific 
ground state, which is needed when using the system as 
memory. 

As a specific example we analyze Ho on Pt(lll) with 
crystal-field Hamiltonian as determined in Ref. [l|, but 
we also compare with the situation obtained for the pa¬ 
rameters of Ref. The quantitative comparison with 
the experiments requires assumptions about a number of 
further parameters, which we explore in the latter part 
of this paper in some detail. For the comparison we also 
present experimental data which resolve, beyond what 
has been reported in Ref. [I|, the lifetimes for each one 
of the low-energy states. 


II. THE MODEL 

A. The crystal-field Hamiltonian and eigenstates 

We investigate the dynamics of a magnetic adatom 
placed on a high-symmetry position on a metal surface 
in a situation where the symmetry stabilizes the degen¬ 
erate ground states against the dominant perturbations. 
Such a situation is realized for Ho on Pt(lll), which 
we consider as a specific example, but our analysis and 
conclusions can easily be generalized to other magnetic 
adatoms. Ho adatoms on a Pt(lll) surface have been in¬ 
vestigated experimentally in Refs. [l|,H,[3- The Ho atom 
has strong spin-orbit coupling, therefore the total angu¬ 
lar momentum is a good quantum number with J = 8, 
leading to 17 states in the multiplet to be studied. A sin¬ 
gle adatom on a crystal structure with trigonal symmetry 
can be described by a crystal-field Hamiltonian adjusted 
to the symmetries [13, Uli , 

Hcf= (1) 

n^2,4,6 

m=0,3,6 

m<n 

= 3 J^ + b2- 35 4 

Here Ojf are the Stevens operators expressed in pow¬ 
ers of Jzi J+ and J_ angular momentum operators of 
the Ho atom [I^. Due to the trigonal symmetry the 
operators J+ and J_ appear in powers of multiples of 
three. Above, the first three terms are shown explic¬ 
itly (with [•,•]+ denoting the anti-commutator), further 
ones are listed in the Appendix. The coefficients 
are crystal-field parameters, which in Ref. [I| were de¬ 
termined from ab-initio calculations and compared with 
experiment. For Ho and many other magnetic adatom 


systems the leading term is oc with a negative coeffi¬ 
cient i?2 < 0- The difference in the conclusions reached 
in Refs. [I| and [3] arise from different values of the pa¬ 
rameter i?4. We first proceed using the values of Ref. [l| 
but we will comment of the situation of Ref. [3] in section 
lYB 

The eigenstates of the crystal-field Hamiltonian Hcf 
can be divided into three families of states, |'0m)i 

and |'0m)- Each one of these states is a superposition of 
different Jz eigenstates with magnetic quantum numbers 
differing by multiples of 3. The lower index m of each 
state, with —8 < m < 8, denotes the dominantly con¬ 
tributing Jz eigenstate. In Fig. [5] we plot the eigenener- 
gies versus the expectation value {Jz)- They are marked 
by circles, squares and triangles for the three families 
-b, —, 0. Since the dominant contribution to the eigenen- 
ergies arises from the first term oc — 4 the energies fie 
approximately on an inverted parabola. Note that by 
plotting the energy versus {Jz) we can also represent ar¬ 
bitrary superposition states of the aforementioned basis 
states, which is useful for the following discussions. 

The two degenerate ground states of the system 
and IV’-s); with differing angular moment pointing into 
or out of the surface, belong to two different families (-|- 
and —). The same applies for the first excited states 
|4) s-iih IV'-r); s-s well as for various higher ones. These 
states are superposition of non-degenerate J^-eigenstates. 
In contrast, the states IV’gs) IV’ea); belonging to the 
0—family, are superpositions made up of degenerate Jz- 
eigenstates, coupled by the operators O4, 0\ and Og. 
They naturally split into symmetric and antisymmetric 
(s,a) combinations, both with vanishing {Jz), and are 
marked accordingly in Fig. [2] a). The same applies for 
14,) and |4„). 



FIG. 2. a) Energies of the eigenstates of the Hamiltonian Hcf 
plotted versus the expectation value (Jz). b) Energy expec¬ 
tation values of the steady states of the quantum master 
equations for strong dissipation due to an applied voltage of 
U = 7.3 mV. The labels refer to the dominantly contributing 
Jz-eigenstates. 
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B. Spin-dependent electron scattering 


From these relations we find ^ since 


For control and read-out of the magnetic state the Ho 
adatom is coupled to a spin-polarized STM tip. The 
electronic reservoirs influence the dynamics of the sys¬ 
tem by the scattering of electrons in the electrodes from 
the atom and by the tunneling of electrons between 
tip and bulk substrate via the atom (see Fig. [T]). Ac¬ 
cordingly the total Hamiltonian consists of three parts, 
H = Hs + Hues + He, describing the system, the elec¬ 
tron reservoirs in the bulk substrate and tip of the STM, 
and the coupling terms, 


Hs 

Hjies 


He 


Hep + gjfJ-sB ■ J, 

i.^ka + eU°‘) elicit, 

a^T,B 



a^OL' —T^B 
kk' 




( 2 ) 

(3) 


(4) 


In ([2]) we added the contribution due to an applied or 
stray magnetic held B, with gj being the Lande-factor, 
gLB the Bohr magneton. The bath electron, with cre¬ 
ation (annihilation) operators c^l (c^^), have energies 
where a can be T or H for the spin-polarized tip or 
the bulk substrate. We also account for the voltage of the 
baths 17“. In (|4|) we concentrate on the interaction of the 
spins of the electrons with the angular momentum J of 
the atom, which leads to scattering and tunneling pro¬ 
cesses, with or without spin flip, with amplitudes 
described by the three terms in He- For a = a' = B 
these terms describe the scattering of a bulk electron, 
while for a ^ a' they describe tunneling between the tip 
and the bulk via the Ho atom. Other couplings of the 
total angular momentum to, e.g., phonons or the radi¬ 
ation field and the resulting dissipative effects are not 
considered explicitly in this work, but will be accounted 
for in a qualitative way in section fV Cl We further have 
to keep in mind that the measured current between tip 
and substrate is also caused by further spin-independent 
coupling terms. 

The eigenstates lipm) of crystal-field Hamiltonian 
have the property that two eigenstates with the same 
energy but from different families have vanishing matrix 
elements (V'mlJi.lV'Im) = 0 for z/ e {+,-,z}, a = ± 
and all m. It is a consequence of two properties of the 
system: First, the symmetry of the adsorption site, 
which leads to the Stevens-operators introduced above 
with operators J+ and J_ appearing in powers of mul¬ 
tiples of three. Second, the time-reversal symmetry for 
B = 0 with the following properties of the time reversal 
operator T (see Ref. Q), 


r^=l, (xk) = (Tz^lTx), 

rj. = -j.r, r\rj = \^z:n)- (5) 


={rrjrj,.^pz:,)* = - 

= - • 

( 6 ) 

Thus, a transition between the two ground states with 
m = ±8 (and similar m = ±7) cannot be induced by a 
single electron scattering, which is crucial for their ob¬ 
served long lifetimes [if. 

On the other hand, at non-zero temperature T yZ 0 
and under the influence of an applied voltage U, scatter¬ 
ing and tunneling of electrons may lead to transitions to 
exited states and eventually to transitions between the 
two ground states. Additionally, time-reversal symmetry 
breaking terms in the Hamiltonian, such as a magnetic 
field, give rise to direct transitions between the ground 
states. These effects lead to a finite relaxation time Ti, 
which has been probed in the experiment and is the first 
and central quantity to be studied in this paper. In a later 
section we will also investigate the decoherence time T 2 , 
i.e., the time scale for the decay of a coherent superposi¬ 
tion of the two ground states. 

The temperature in the experiments [I| was as low as 
0.7K (i.e., ksT Ri 0.060meV), whereas the typical sys¬ 
tem energy, i.e. the first excitation energy is 7.7 meV. 
Thus, the system is clearly in the quantum regime. 


III. QUANTUM MASTER EQUATION 

A. Damping due to electron scattering and 
tunneling 


The description of the system in the quantum regime 
under the influence of the electronic reservoirs requires 
solving the reduced quantum master equation appropri¬ 
ate for open quantum systems [H, [l^ ; 

p{t) =i[p{t),Hs]+ [ dt'i:{t-t')p{t'). (7) 

Jto 

Here and below we set h = 1. The influence of the two 
electronic reservoirs enters in the dissipative kernel S(t— 
t'). We assume it to be of a Lindblad form and use a 
Born-Markov approximation. In the interaction picture 
the quantum master equation for the 17 states of the 
system then reduces to 

poo 

pi= - ^ / dt' (8) 

/ , Jo 
=+, — ,2 

a,a' =^T,B 

I [Ju{t)J,,'it')pi{t) - Ju'{t')pi{t)J,,{t)] (7““, (< - t') 

+ [piit)j,.:{t')Mt) - Mt)pi{t)j,,{t')] cz'if -1) I, 
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with the dissipative kernel expressed by the correlation 
functions 


CZ' it) = E \"i^tkUt)st''k.' (0)) (9) 

k,k‘ 


with 


„aa _ „a' aa' _ at„a' 

®fcfc'- ^kk'+ ‘-fc-f'fe't’ 

^kk'z — Cfc-|-Cfc/^ - C^4,Cfe/^. 


( 10 ) 


Assuming « t““ and introducing the spin- 

dependent electron densities of states at the Fermi-edge 
iV“ with a =tii we obtain the Fourier transformed of 
the correlation functions, 


C““'(A„^) = |t““'|2Ar“ArAc(A„„ + e(C/“ - C/“')) 
Cr+iK^) = |r“'|2iV^“iVf'C(A„„ + e(C/“ - C/“')) 
C““'(Anm) = |t““? (iVf iVf' + 

xC(A„™ + e([/“-[/“')). (11) 


They are evaluated at the energy differences of the atomic 
system, A„m = — E^ shifted by the applied voltages. 

Here we introduced 

C(„) = / /(E) [1 - /(E - „)] diJ = 

where /(E) = [e^AfcsT) _|_ ^j-i jg Fermi function. 

As usual in the context of the tunneling magneto¬ 
resistance we define the tip polarization rj = {P^ — 
Ei.)/iEt + -Pj.) G [~li 1]; where the spin up/down popu¬ 
lations are proportional to the densities of states oc 
= N'^ ■ ^(1 ± 77 ). The bulk electrode is assumed to 
be non-polarized, hence . The remaining pa¬ 

rameters, apart from the polarization 77 , can be lumped 
in the coefficients 

In the following discussions of the relaxation pro¬ 
cesses we will concentrate mostly on applied voltages 
U = [/^ — of the order of or exceeding 3 mV. In 
this case the effect of tunneling electrons is stronger than 
that of scattering electrons in the bulk electrode, which 
in turn is assumed to be stronger than the scattering in 
the tip. We therefore set ctt = 0 in most of the paper. 
Furthermore, we start with cbb = 0, but we will ana¬ 
lyze the influence of bulk electrons scattering in section 
IV Al Also for the calculation of the dephasing time T 2 for 
zero current in section ED the scattering terms need to 
be taken into account. 


B. The proper basis states and 
environment-induced superselection 

If the system behaves sufficiently classically, it is possi¬ 
ble to reduce the quantum master equation to rate equa- 



FIG. 3. Expectation value (Jz) versus the voltage U for a tip 
polarization r] = —0.15, and in the inset versus the polariza¬ 
tion 77 for U = 7.3 mV. The other parameters are Bz = 10“® T 
and ctb = 3.41 ■ 10®(meV s)“^. 


tions for the populations Pi of the different states 
d P 

^ = E i^pPj - , (12) 

3 

with Tji describing the transition rates between the dif¬ 
ferent states. Such a reduction is possible if the coher¬ 
ences, i.e., the off-diagonal components of the density 
matrix decay much faster than non-equilibrium popula¬ 
tions. The question remains, what is the appropriate 
basis to decouple the coherences and populations. One 
frequently used option for this basis are the eigenstates of 
EIcf- It is the appropriate choice if the influence of the 
baths is weak, and the eigenstates are only little affected 
by their presence. This choice and the resulting reduc¬ 
tion to rate equations has been used, e.g., in Refs. [IMl 
and allows an easy interpretation of transition channels. 

In general, the quantum master equation ([S]) for the 
density matrix can be rewritten in the form p(t) = Mp(t) 
where the density matrix p (here 17 x 17) is arranged as 
a vector p (with 289 components) and all the dynamics, 
coherent and dissipative, are included in the components 
of the large (289 x 289) matrix M (for details see Ap¬ 
pendix). The spectral decomposition of this matrix, 

A4 Pn — kJln Prif (1^) 

gives access to various relaxation rates. It also yields the 
steady-state populations p{t —>■ oo) = pst = Po, which 
is the eigenvector corresponding to the zero eigenvalue 

TTlo = 0 [2^. 

If the effect of the baths is not weak a reduction to a 
classical rate equation is not allowed. Instead we should 
study the eigenvalue problem of Eq. (USD. Its eigenvec¬ 
tors are superpositions of different eigenstates of EIcf- 
This is in particular true for the 0-family. We recall 
that due to degeneracies the eigenstates of EIcf corre- 
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spending, e.g., to m = ±6 are the symmetric and an¬ 
tisymmetric coherent superpositions of (predominantly) 
the two Jz basis states, resulting in a vanishing expecta¬ 
tion value (Jz) displayed in Fig. [5] a). Due to the dissi¬ 
pation this coherence may get destroyed. This scenario 
has been termed “environment-induced superselection” 
principle by Zurek [3, 0. Decoherence selects favored 
“pointer” states which are stable under the influence of 
the environment [2l|. These states are the steady-state 
eigenstates of the full quantum master equation which we 
will denote as \(j)m) with m = —J, —J -|- 1,..., J. If the 
dissipation is strong these eigenstates are actually much 
closer to the original Jz basis states with non-vanishing 
expectation values (Jz) ~ ±6. Similarly the states cor¬ 
responding to m = ±3 get modified. The situation with 
strong dissipation is depicted in Fig. [2]b). 

It is instructive to study the transition between the 
two limiting cases. We focus on the two stationary states 
of the doublet \(jy+e) and |(/>-6) and calculate their Jz ex¬ 
pectation value. The result is depicted in Fig. [3] For 
weak dissipation, e.g., for small U the states reduce ap¬ 
proximately to the Hcf eigenstates |'0g^) and IV'ea) 

( Jz) ~ 0, while for strong dissipation they approach the 
states with maximum expectation values (Jz) ~ ±6. For 
low voltages the electrons do not have enough energy 
to scatter into higher excited states and the superposi¬ 
tion states remains stable. For higher voltages exceeding 
~ 6 mV and nonzero temperature excitations are possible 
and the superposition is destroyed. In Fig. [3] we further 
note a special feature, namely a dip at U = 7.7 mV, since 
at this voltage the electrons are in resonance with the first 
excitation energy. 

In the inset of Fig. [3] the dependence of the 
environment-induced superselection on the polarization 
of the tip is shown. Without polarization, the states have 
no preferred basis in which they evolve upon scattering 
of the electrons. But already a very small polarization in 
z-direction (mind the scale of the axes) is enough to drive 
the system to the pointer states with maximum (Jz). 

The environment-induced superselection strongly in¬ 
fluences the relaxation processes, which will be studied 
in the following section. 


IV. RELAXATION TIME Ti 

If the system is mainly in one of the two degenerate 
ground states, the relaxation time Ti towards the steady 
state is given by the smallest non-zero eigenvalue of 
the matrix A4 (see Eq. m), corresponding to the 
eigenvector pi = (1,0,...., 0, —I)^. (Here we assumed 
an ordering such that the hrst and last entries of the 
eigenvector are the ground states). The inverse of the Ti 
time, 1/Ti = —mi = r_8_j.+g-|-r+g_>_g, accounts for all 
relaxation channels from one ground state to the other. 
Note that the switching rate F_g_i.+g from |())_8) to 
\(j)8) accounts for the direct transition but also all those 
via excited states. It thus differs from the rate r_gg of 


the rate equations (HD), which describes only the direct 
transition. The relation between the two switching rates 
F_g_>+g and r+8_)._g (and related lifetimes) follows from 
the steady state populations P-s and P+g of the two 
ground states, r.g/r+s = r+8_j._g/F_8^+8 = P.g/P+g. 


A. Voltage dependence 

To compare with the experiments of Ref. [I[ we need 
to determine the coupling strength for tunneling of the 
reservoir electrons via the Ho atom. For this purpose we 
calculate the current Ith from the dissipative part of the 
master equation (for details see Appendix). By compar¬ 
ing with the experiment we should be able to determine 
the coupling strength ctb- However, the current Ith de¬ 
scribes only the current where electrons scatter due to the 
spin - angular momentum interaction. The total current, 
which in the experiments was always kept at I Exp = 1 nA, 
includes a ‘leakage’ current I Exp = Irh + I Leak- It 
can arise due to electrons tunneling directly between tip 
and bulk or due to electrons scattering with shells other 
than the 4/ shell which forms the basis of the consid¬ 
ered angular momentum states. By comparing the cur¬ 
rent which involves spin flips, and accordingly depends 
on the the spin state of the Ho atom, with the state- 
independent current we get a rough estimate. For the 
following discussion we assume that Ith amounts for 
roughly 10% of the total current Iexp, i-e. Irh = 0.1 nA. 
For U = 3 mV this is obtained for a coupling strength of 
Ctb = 3.41 • 10®(meV s)“^. 

We are now ready to analyze the dependence of the 
relaxation time on the applied voltage. In Fig. |4] a) we 
plot the results for the temperature T = 0.7 K chosen 
in the experiments, and for comparison also at several 
higher temperatures T = 1.4K, T = 3.5K, and T = 7K. 
We compare the Ti time, as obtained from the numerical 
solution of the full quantum master equation, and the 
result obtained from the approximate rate equations m 
in the basis of Hcf eigenstates. Similarly we compare in 
Figs, mb) the Jz expectation values of the steady state 
solutions of the master equation and those of the Hqf 
eigenstates for two different values of U. For reasons of 
numerical stability we include in all our simulations the 
effect of a very weak magnetic field applied in z-direction 
(Hz = 10“®T). Otherwise the two ground states get 
completely decoupled within our numerical precision, and 
divergences appear, or the reduced density matrix is no 
longer positive semidefinite. 

Focusing on low temperature, T = 0.7K, we note 
that for low voltages, U < 3 mV, the six lowest lying 
steady states of the full quantum master equation have 
very similar properties as the Hqf eigenstates. In this 
regime transitions between the two ground states j'l/'^) 
to |'0Zg) are caused mostly by the (weak) symmetry¬ 
breaking magnetic field and are thus voltage indepen¬ 
dent. The full master equation yields shorter lifetimes 
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FIG. 4. a) Relaxation time Ti and decoherence time T 2 of 
the two ground states versus the applied voltage U obtained 
from the full quantum master equation (1131) (full lines) and the 
rate equation m based on Hcf eigenstates (dotted lines). A 
weak magnetic field Bz = 10“® T is applied, the temperatures 
are T — 0.7K, T = 1.4K, T — 3.5K and T = 7K, and the 
coupling strength is ctb = 3.41 ■ 10®(meV s)“^. b) Energy 
expectation values versus the expectation values {Jz) for the 
steady states of the full quantum master equation and for the 
Hcf eigenstates at T — 0.7 K for U = 3 mV and U = 7.3 mV. 
c) Relaxation time Ti versus the coupling strength ctb as 
obtained from the full quantum master equation and the rate 
equation for T = 0.7 K and U — 7.3 mV. 


than the rate equation in this regime, because of addi¬ 
tional coherent transitions contained in the full theory. 

At higher voltages we observe for T = 0.7 K and 
T = 1.4 K in (the semi-log plot of) FigUla) an exponen¬ 
tially activated behavior. In this regime the full master 
equation and the rate equation yield very similar results. 
Here the main switching channel is via the first excited 
states, 107 ) or | 0 - 7 ), followed by a fast decay to the other 
ground state on the other side of the parabola. Since 


the first step is the bottleneck of the process we have 
l/Ti Ri Egy -I- r_ 8 _ 7 . For an estimate we ignore the ef¬ 
fect of the tip polarization and of a magnetic field (i.e. 
Fgr r; r_8-7) and get 

r87 =ctb\ P • C{E7 -Es- eU) (14) 

(Er-Es-eU) 

~ 16 cts ^(^Ej-Es-eU)/kBT _ I 

For still higher voltages, U > 7 mV, the results ob¬ 
tained in the two approaches differ again. In this regime, 
in the frame of the rate equation excitations to the sym¬ 
metric and anti-symmetric i/cj’-eigenstates |'0g^) and 
|0g£,) become possible. These two states have high tran¬ 
sition rates between each other and thus provide a ‘short¬ 
cut’ for the decay. However, in the full master equation 
these superpositions states are replaced by the steady 
states 10+6) and |0-6), with properties illustrated in Figs. 
[ 3 ] and |4] b), which in this parameter regime are actually 
close to {Jz) = ±6 states. They are weakly coupled, 
and the shortcut is no longer open, which increases the 
relaxation time. 

For even higher voltages, U > 8.6 mV, tunneling elec¬ 
trons have enough energy to overcome also the second 
energy excitation gap of AEq^ = Eq — E^j k. 8.6 meV, 
which is the largest gap of the system. From this point 
on, sequential scattering over the top of the parabola is 
the main transition channel, and the slope of Ti(U) ver¬ 
sus U changes. 

At higher temperatures, the different regimes get 
smeared out, as can be seen in the plot of Fig. SI 
especially for T = 3.5K and T = 7K. In this regime 
the main transition channel is always via higher excited 
states. For voltages larger than U r; 9 mV the results 
for all temperatures are very similar, because in all cases 
most of the electron scattering leads to the transition 
over the top of the parabola. Again, the rate equations 
overestimate the role of transitions via the ’shortcut’ 
states 10-g^), |0g„), |0§^) and |0§„). 

In Fig. |3]c) we investigate the dependence of the life¬ 
time on the coupling strength ctb for low temperature 
and U = 7.3 mV. From the rate equations we find sim¬ 
ply that Ti time decreases proportional to IjcTB- But 
the solution of the full quantum master equation yields 
different results. In the considered regime excitations 
to the 10 + 5 ) and |0-6) play a role. The coherence lead¬ 
ing to these superposition states is increasingly destroyed 
with growing coupling strength. As a result in the range 
10®(meV s)“^ < Ctb < 10^(meV s)“^ the stronger cou¬ 
pling even stabilizes the system by decoupling the two 
sides of the parabola. 


B. Further details of the experiments 

To compare with the experiments of Ref. [I|, two 
further effects have to be taken into account. First, 
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FIG. 5. Relaxation time Ti and decoherence time T 2 of the two ground states versus the applied voltage f/ in a weak magnetic 
field Bz ~ 10“® T at T = 0.7 K for the full master equation solution and the rate equation m with the Hcf eigenstates, 
a) With and without modulation voltage Umod = 0.8 mV. b) With and without tip distance correction. For comparison the 
experimental data are shown in red. The error bars are too small to be visible. 


in the experiment the applied voltage was modulated, 
Utotit) = U + V2UmodCOs{ujt), with Umod = 0.8mV. 
The frequency of the modulation, w = 720 Hz, is several 
orders of magnitude faster than the timescales of inter¬ 
est. We model this fast oscillation by suitably averaging 
the correlation functions with the distribution h{U') = 
l/^[ 2 (eC/„od) 2 -(eCp 2 ]-i /2 |^,| < ^Umod and 

h{U') = 0 otherwise [^. This means, e.g. 

C^-{^nm) = + v) 

X J dU'CiAnm + eU^+ eU')h{U'). (16) 

The effect on the lifetimes is shown in Fig. [5] a). In 
essence, the modulation amounts to a shift U ^ U + 
■\/2 ■ Umod, by approximately 1.1 mV. 

Second, in the experiment the voltage-induced tunnel 
current was measured, and in fact, by adjusting the dis¬ 
tances between tip and electrode, was kept fixed at the 
value I Exp = 1 nA. To keep the current fixed for differ¬ 
ent voltages U, we allow in the simulations ctb{U) to 
be voltage dependent. We checked numerically that the 
tunnel contact behaves close to Ohmic, Ith ~ U/Rtuu- 
We therefore adjust the coupling ctb{U) oc 1/C/ to keep 
the current constant. We recall that an estimated 90% of 
the current is leakage current. With ctb{U = 3mV) = 
3.41 • 10®(meV s)“^ we arrive at Irn = O.lnA. We will 
proceed using these values in all simulations reported be¬ 
low. All this being said and done, we note that the effect 
of the adjustment, which is included in Fig.[S]b), is weak. 


V. DEVIATIONS FROM THE IDEAL 
SITUATION 

When comparing the calculated Ti times with the ex¬ 
perimental data, as illustrated in Figs. [5l we note that 
the theory produces far too long times. Therefore, we 
need to take a closer look at the experiment and possible 
deviations from the ideal situation assumed so far. 


A. Scattering of bulk electrons 

Up to now we ignored the effect of bulk electrons scat¬ 
tering from the Ho atom, and the question arises whether 
it could be the source of the mismatch between the¬ 
ory and experiment. The scattering processes are eas¬ 
ily included in the quantum master equation, and their 
effects are illustrated in Fig. [51 In the first panel we 
show how the relaxation time depends on the coupling 
strength cbb for a fixed value of the tunneling strength 
ctb = 3.41 • 10®(meV s)“^ and voltage U = 5 mV. 
As long as cbb is smaller than ctb the lifetime re¬ 
mains nearly unchanged. For stronger cbb, the scat¬ 
tering of bulk electrons leads to suppression of coher¬ 
ent transitions and thus to longer lifetimes. This arises 
because of the low temperature of the bulk electrodes, 
which cools the system into the ground states, whereas 
the tunneling electrons due to the applied voltage have 
enough energy to excite the system. The combination 
of voltage-dependent tunneling and voltage independent 
scattering is illustrated in Fig. [5] b) (including the ef¬ 
fects of modulation voltage broadening and tip distance 
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FIG. 6 . a) Relaxation time Ti versus the bulk electron scattering strength cbb as obtained from the full master equation m 
and the rate equation m based on Hcf eigenstates. The parameters are U = 5 mV, Bz = 10 ® T, ctb = 3.41 • 10®(meV s) 
b) Relaxation time Ti versus the voltage U for the model with modulation voltage broadening and tip distance correction with 
and without bulk electron scattering strength cbb ~ 1-0 ■ 10^ (meV s)“^. 




FIG. 7. Relaxation time Ti of the two ground states under the influence of deviations from the perfect situation. We compare 
the results obtained from the full master equation (O and the rate equation m based on the Hcf eigenstates, a) Ti versus 
Stevens parameter B 2 characterizing the deviation from the perfect trigonal symmetry, b) Ti versus the strength of Gaussian 
broadening 7 chosen to account for further lifetime broadening effects. In both plots we choose the parameters U = 5 mV and 
Bz = 10~® T. The insets show the Jz expectation value of the steady states and the Hcf eigenstates. 


correction described above). We assume a bulk elec¬ 
tron scattering strength cbb = 1-0 • 10^ (meV s)“^ which 
is slightly higher than the tunneling coupling strength 
Ctb = 3.41 • 10®(meV s)“^. Again we note the increase 
of the lifetime as a consequence of the scattering. 

B. Breaking the Gs^-symmetry 

As an example of a symmetry breaking term, we con¬ 
sider the effect of the Ca^-symmetry-breaking term of 
the Stevens operators, B^-O^ = B^- [Jz, J+ + J-]+ [^ . 
This term arises if the tip is not perfectly centered over 
the Ho atom, or if nearby surfaces or imperfections in 
the crystal break the symmetry. It breaks all the rota¬ 


tional symmetries of the system. A magnetic field in the 
xy-plane would have a similar effect. 

The symmetry-breaking parameter B 2 is varied in 
Fig. [7] a) between 10“^°meV and 10“^meV, which is 
still orders of magnitude lower than the leading crystal- 
field term = —0.239 meV. As a result of the sym¬ 
metry breaking the eigenstates of Hcf get mixed, and 
the protection against direct transitions is lost. In Fig. 
[7] a), the voltage is chosen to be C/ = 5meV, so the lead¬ 
ing transition is directly between the two ground states. 
With rising strength of the relaxation time decreases 
Ti oc {BI)~^. In the frame of the rate equations the 
switching arises since electron tunneling directly couples 
the two states. In the frame of the full master equation 
scattering of electrons destroys the phase coherence of 
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these superpositions and the resulting expectation values 
(J^) of the steady states are closer to the parabola, char¬ 
acteristic for the Hcf eigenbasis (see the ground states 
of the inset of Fig. Ha)). This means, that the symmetry 
protection gets restored and the lifetimes increase dras¬ 
tically. All this depends on the coupling strength ctb, 
which is the important parameter for the superselection 
(see Fig. He)). Unfortunately, the parameter is not 
independently accessible in the experiment, and it is dif¬ 
ficult to draw more precise conclusions. 

C. Noise in the circuit 

To account for the influence of further perturbations, 
such as thermal noise in the electronics, we introduce 
a lifetime broadening. Specifically we average the re¬ 
sults obtained so far, assuming a Gaussian broadening 
of the tunneling electron energy distribution g{SE,^) = 
Xj exp[—( 5 U^/( 27 ^)] with a width characterized 

by the parameter 7 [23| . Its influence on the lifetime T\ is 
visualized in Fig.Hb). The relaxation time is strongly re¬ 
duced by this broadening because the broadening allows 
excitations into one of the or \(j)-^) states, followed 
by a subsequent decay to the other ground state. 

By fitting the parameters characterizing the two devia¬ 
tions from the ideal situation we manage to obtain results 
for the Ti times close to the the experimental ones (see 
Fig. [HI). We plot results for two values of the tempera¬ 
ture T = 0.7 K and T = 1.4 K, where the first one is the 
bath temperature in the experiment, while the second is 
chosen to account for some electron heating induced by 
the current. (An increased temperature alone would not 
be sufficient to explain the discrepancy between the sim¬ 
ulations and the experiment.) Although with the fitting 
we reached a reasonable agreement with experiment we 
have to acknowledge that the fit is not conclusive, since 
the number of data points is too low to determine the 
parameters independently. 

D. Magnetic field dependence 

Next, we investigate the effect of an applied or stray 
magnetic field which is probably present in the exper¬ 
iments. Its influence depends strongly on the values of 
the other parameters. In Fig. [9] a) we show the resulting 
modification for the regime where the main transition, 
although with small rate, is directly between the states 
|'0^) and IV’Ig)- choose U = 5 mV, i.e., we are still in 
the regime of volt age-independent relaxation time of Fig. 
Ha). As shown in Fig. Ha) the lifetime is strongly re¬ 
duced by the symmetry breaking magnetic field , since 
the symmetry protection of the two ground states is lost. 
This behavior is obtained both from the rate equation 
and the solution of the quantum master equation. 

Interestingly, as illustrated by Fig. |9|b), the behavior 
can be completely different for different parameters. In 



FIG. 8 . Relaxation time Ti and decoherence time T 2 of the 
two ground states versus the applied voltage U at T — 0.7 K 
and T = 1.4 K for the modulation amplitude U-mod ~ 0.8 mV, 
lifetime broadening 7 = 0.95 meV, magnetic field = 
1 ■ 10“® T, symmetry breaking B 2 = X ■ 10“"^ meV, and tip 
spin polarization rj — 0.15. We compare the results obtained 
from the full master equation (1131) (solid lines) and the rate 
equation (HD based on the Hcf eigenstates (dashed lines). 
The theoretical values are compared with the state-dependent 
data from the experiment with error bars indicating the sta¬ 
tistical errors of the measurement. 

this example the lifetime increases when a field is applied. 
For the chosen parameters the B^ field stabilizes the Jz 
eigenstates and reduces the switching, which otherwise 
would be induced by the symmetry breaking term B^. 
We conclude that a detailed study of the magnetic field 
dependence should provide a better understanding of the 
different perturbations acting on the Ho adatom. 


E. Alternative choice of the crystal-field 
parameters 

In a recent publication Donati et al. Q reported about 
x-ray absorption spectroscopy (XAS) and magnetic cir¬ 
cular dichroism (XMCD) measurements performed with 
Ho on Ft (111). The results could be fitted with a crystal- 
field Hamiltonian of the form m but with only the co¬ 
efficients B 2 = — 140^eV and B^ = 1 ^eV differing from 
zerm which differs significantly from what Miyamachi et 
al. []| used, and what we assumed so far in this paper. 
To simulate the lifetimes in our theory with these pa¬ 
rameters, we add a small Bf « 0.3 ^eV which we take 
from Miyamachi et al. [I| , because otherwise the ground 
states decouple completely. With these inputs, we ar¬ 
rive at the level scheme shown in the inset of Fig. [101 
Most important we note that the ground states are now 
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FIG. 9. Relaxation time Ti of the two ground states versus a magnetic field Bz as obtained from the full master equation m 
and the rate equations m based on Hcf eigenstates. The insets show the expectation value of Jz for the two descriptions, 
a) Results for the ideal model for U = 5 mV. b) Results with symmetry-breaking and Gaussian broadening corresponding to 
the fits of Fig. [ 8 l i.e., Umod = 0.8 mV, 7 = 0.95 meV, and B 2 — 4: ■ lO”'^ meV. 



FIG. 10. Relaxation time Ti of the two ground states versus 
the voltage U as obtained from the full master equation and 
from the rate equations m based on Hcf eigenstates for 
Bz = 10“® T. The inset shows the expectation value of Jz for 
the two descriptions for U = 8 mV. The black line shows the 
function f{Jz) = — 140/ieV • O 2 + 1 fJ-eV ■ O4 -|-const.. 


the IV'es) IV' 6 a)i which are much stronger coupled 
(not symmetry-protected) than the states and \ipZs) 
and hence should have much shorter lifetimes. We ana¬ 
lyzed the relaxation rate in the same way as for the other 
model, with results shown in Fig. [TUI We note that the 
solution of the rate equation yields very short relaxation 
times of the order of nanoseconds. On the other hand, the 
simulation based on the full master equation yields longer 
lifetimes. The reason is again the environment-induced 
superselection which destroys the superpositions. The 
differences between the two sets of eigenstates are pro¬ 
nounced, as can be seen in the inset of Fig. [TUI hi fact the 
difference between the two theoretical approaches is even 


more pronounced than found in the model based on the 
parameters of Miyamachi et al. We further note from 
Fig. [TU| that higher voltages even stabilize the ground 
states, because the excited states couple less to the states 
on the opposite side of the parabola than I'l/'+e) and |' 0 ° 5 ) 
do, and hence the excitation of those states reduces the 
transition rate. This voltage dependence is in stark con¬ 
trast to the observations made by Miyamachi et al. [I|, 
where the lifetimes decrease with increasing voltages. 

The comparison of theory and experiment on the 
voltage dependence of the relaxation time supports the 
Stevens parameters used by Miyamachi et al. [J. They 
were derived by ab initio DFT calculations for a situa¬ 
tion where individual Ho atoms were adsorbed on high- 
symmetry fee sites on the surface of Pt(lll). In contrast, 
the experiments of Donati et al. @ were performed with 
a high coverage of Ho atoms of 0.04 monolayers occupy¬ 
ing a mixture of hep and fee sites. Further investigations 
are needed to clarify whether this is the origin of the 
differing results. 


VI. DECOHERENCE TIME T2 

The decoherence time T 2 is the time scale on which 
the phase information in a coherent superposition, here 
specifically of the states and I'i/’Zg), gets lost. In the 
quantum master equation treatment of the problem, T 2 
is obtained from the decay rate of the off-diagonal ma¬ 
trix element p+g-g, which is given by the corresponding 
matrix element in the matrix M, 

I/T 2 = — A4g_g_>g_g 

~4cts {'tptlJzliJs) {'^ZslJzliJZs) C{-eU) 

Ri 264 CTseH. (17) 




































11 


The resulting T 2 times are plotted for different parame¬ 
ters in Figs, m [5] and HI It turns out that the T 2 time 
depends mostly on the current, i.e., on the number of 
scattered electrons. Roughly one can argue that each 
scattered electron dephases the superposition state, in¬ 
dependent of the energy of the electron. Thus, we detect 
only in Fig. |3]a) and Fig. [S]a) a dependence of T 2 on the 
voltage U, because in all other plots the current is kept 
constant (achieved by the adjustment of ctb{U)). We 
note that the decoherence times T 2 are always very short 
(~ 10 “^*^ s) making the considered memory unsuitable 
as a quantum mechanical bit (qubit). Additionally, we 
found that T 2 does not depend on the symmetry breaking 
B2, the broadening 7 , or the magnetic field B^- If the 
voltage is set to zero, all terms (±A„m) are similarly 
important. Hence T 2 is given by 

1/22 = — Afg-S-s-S-S 

~4(^ c„a') C(0) 

acn' 

m 264 kgT E (18) 

Ota.' 

If one assumes that J2aa' ^aa' ~ 10®(meV s)“^ and the 
temperature is T = 1 K, the decoherence time is roughly 
10 -®s. 


VII. INITIALIZATION 

We have seen that the relaxation time depends strongly 
on various parameters. In this section we will demon¬ 
strate that by switching parameters suitably we can ini¬ 
tialize a specific angular momentum state |24| . This en¬ 
ables the writing process if the Ho atom is used as a mem¬ 
ory. As a specific example we study what happens to the 
state IV'Zg) when we pulse the system for a time tp with a 
voltage U and then let the system relax for 1 fis without 
a voltage applied. After the relaxation, the population 
of the state \tps) is measured, which provides the infor¬ 
mation about the switching probability S-8^8{U,tp). 

In Fig. [TTl the switching probability S-8^8{U,tp) as 
obtained from the quantum master equation is plotted 
versus the strength of the applied voltage pulse for differ¬ 
ent spin polarizations 77 of the tip. (We neglect again the 
influence of the bulk electron scattering.) For U > 20 mV 
the value of S-8^8{U,tp) may get close to 1 , i.e., the 
pulse flips the Ho atom with a high probability into the 
state IV'^). On the other hand, a pulse with opposite 
voltage U < —20 mV leaves the state with high prob¬ 
ability in the state IV'Ig)- Thus, both ground states 
can be prepared by applying voltage pulses with either 
sign. For voltages between —10 mV ^ U < 10 mV ba¬ 
sically no switching of the Ho spin state is induced, 
S-8^8{U,tp) Ri 0. Around U = 20mV the switching 
is strongest as long as rj ^ 0. For an unpolarized tip, 
ry = 0 , a pulse with high voltage produces a balanced 
population of the two ground states. 


The upper inset of Fig. [TTl shows the dependence of the 
switching probability S-8^8{U,tp) on the spin polariza¬ 
tion rj for optimum conditions 17 = 20 mV and pulse time 
tp = 2.5 • 10“^ s. As one could expect, the higher the po¬ 
larization, the better the initialization. But, the switch¬ 
ing between the ground states also requires a certain 
number of electrons, which increases with longer pulse 
time tp, as displayed in the lower inset of Fig. [TTJ For 
the optimal voltage of [/ = 20 mV we calculate a current 
of Irh ~ 0.75 nA. A pulse time of tp r; 250 ns implies 
then that around 1100 electrons are needed to prepare 
the atom in one state. By changing the tip distance in 
the experiment, this dependence could be probed. 

Our simulations suggest that it is possible to prepare 
the system in a required state using pulse times of a few 
hundreds nanoseconds with high fidelity as long as the 
tip polarization is large enough. The required pulse time 
depends strongly on the parameter ctb- Thus, an experi¬ 
ment with different pulse lengths could help to determine 
this value. 



U(mV] 


FIG. 11. The probability for switching from I'i/’Zg) to IV'g") 
versus the applied pulse voltage for different values of the tip 
polarization rj and a pulse time of tp — 2.5 • 10“^ s. The 
upper inset shows the switching probability versus the tip 
polarization for pulse strength U = 20 mV and pulse time 
tp = 2.5 • 10“^ s. The lower inset shows the dependence on 
the pulse time tp for pulse strength of 17 = 20 mV and 7 = 0.2. 


VIII. CONCLUSION 

When investigating the dynamics of the angular mo¬ 
mentum states of Ho on Pt(lll) we found that in an im¬ 
portant regime of parameters the system behaves deeply 
quantum mechanically and cannot be described by a rate 
equation for transitions between the eigenstates of the 
crystal-field Hamiltonian Hcf- Rather the system has 
to be treated by the full quantum master equation. Its 
steady-state basis in general differs from the Hcf eigen¬ 
states, which is an example of the “environment-induced 
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superselection” principle [ 3,0 ■ We analyzed how the re¬ 
laxation time Ti depends on various parameters of the 
system. In ideal situations at low temperatures it should 
be extremely long. In contrast the decoherence time T 2 
is always very short. We further described a method to 
initialize the system in one of the two ground states by 
suitable voltage pulses. 

By a proper choice of the parameters describing devi¬ 
ations from the ideal situation we could roughly fit the 
experiments, with still remarkably long lifetimes. Our 
analysis shows that if the parameters could be improved 
and controlled more precisely the system would acquire 
even longer lifetimes than observed already. The system 
therefore promises to be useful as a single-atom memory 
with the possibility to write into the memory by very 
short pulses of electric currents with high fidelity. 

At this stage, we have to conclude that there are too 
few experimental data available to determine the param¬ 
eters independently. A further detailed investigation of 
the coupling parameter ctb and the other parameters is 
required to identify the main perturbation which limits 
the lifetime of the Ho adatom in the experiment. The 
lifetime depends strongly on an applied magnetic field. 
The analysis of this effect would help obtaining the 
missing information on the parameters. 
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Appendix A: Stevens operators and parameters 
For convenience we list here the Stevens operators 


which are needed to describe the system Ho on Pt(III) 

[ii,[3 

0°=3J2-J(J+1), (Al) 

=35 - 30 J( J -f 1) J2 -h 25 - 6 J( J -f 1) 

+ 3J^{J + 1)\ (A2) 

Ol=^[MJl + Jl) + {Jl + Jl)J,], (A3) 

=231J® - 315 J(J -k 1)4 
-k 7354 + 105J2(J -k - 525J(J -k 1)4 

-k2944-5j3(J-kl)^-k 
40j2(J-kl)2-60J(J-kl), (A4) 


Ol [(114 - 3J( J + 1) J, - 59J,)(4 + 4) 

+ (4 + 4)(ll4 - 3J(J + 1) J. - 59J,)] , 

(A5) 


Ol=\[jl + Jl]. (A6) 

We also list the parameters which were obtained from 
ab-initio simulations and listed in Ref. [l|. 


Anisotropy constant 

value 


-239 peV 

Bl 

86 neV 

Bl 

293 neV 

Bl 

0.186 neV 

Bl 

-1.967 neV 

Bl 

0.630 neV 


TABLE I. Anisotropy parameters as used in Ref. Q. 


Appendix B: Setting np the matrix A4 

In this section we present the approach we use to solve 
the quantum master equation in Born-Markov approxi¬ 
mation. By using the explicit time dependence of the sys¬ 
tem operators Jv(t) in the interaction picture we evaluate 
the time integrals inclu ding the correlation functions of 
the master equation [H, [23 

poo 

^^(±4= / (Bl) 

Jo 

In the eigenbasis oi Hs \n) = En |n), the matrix elements 
of 5'““, (±t) become 


pOO 

{n\ ±tau) \m) = drC'“4(±4 

^0 

^-iHst gZffsT 1^^ ^g2) 


41 e” 


= {n\ Jj./ |m) / dr (7““, (±T)e 

Jo 


iAnmT 


(B3) 


= (n| \m) 


2 J 27r Anm ^ ^ 


(B4) 


Here, Anm = Em — En are the energy differences and 
the tilde over the correlation functions indicate Fourier 
transforms. The imaginary parts of this equation vanish 
in the Born-Markov approximation because in the master 
equation complex conjugate terms are summed. 

For a formulation that can be implemented efficiently 
in source code, the quantum master equation is rewritten 
in the form 


J^P = Mp. (B5) 

Here p = vec{p} denotes the column-vectorization of the 
matrix p, meaning that the (i-kl)*^*' column of the matrix 
is written below the other i ones. The supermatrix M 
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has the dimension dim(p)^, which implies an extension 
of the complexity, but on the other hand the solution of 
this equation is possible with standard numerical tools. 
With the help of the relation 

vec{AXB} = (A 0 B^)vec{X}, (B6) 

where A, X and B are matrices, the transformation of 
quantum master equation into the wanted type is pos¬ 
sible. The symbol 0 is the Kronecker-product of the 
matrices is defined as 


A®B = 


axiB ■ ■ ■ ainB 


(B7) 


JlmlB * ■ ■ 

For the problem considered in this paper we get 


Xi = A4c + Xt D 
Me = i (1 0 FTJ - i7s ® 1) 


Xin = — 


E 


(j.5r/(+r) 0 1 


a,<y'—T,B 


(B8) 

(B9) 


- 0 (J.f) 

+ (l 0 - (J. 0 [SZ'(-T)f)' , 

(BIO) 


where Me contains the coherent part of the master equa¬ 
tion and Md the dissipative part. 


Appendix C: Calculating the Current 

The current is given by the time-derivative of the num¬ 
ber of particles of the tip Nxit) = Yhka 

lTh{t) = = ie{[H,NT{t)]). (Cl) 

In the commutator [iJ, iVj’(t)] only the coupling Hamil¬ 
tonian He survives, which leads to 

iTh = {{J+c^iclA + 

kk' 


{Jz 


Bf T Bf T 


) 


l J [r^t„B _ Tt B 


(C2) 


This form of the current bears similarity to terms in the 
master equation. It contains all the tunneling events from 
the tip to the bulk with a positive sign and those in op¬ 
posite direction with a negative sign. Proceeding similar 
to the steps used for the dissipative part Md oi the mas¬ 
ter equation in Born-Markov approximation^ and concen¬ 
trating on the stationary limit one finds [28l | 




-fr/i =Tr IthP = e'^{MA)(i-i)-^'^+oAPSt)i- 


(C3) 

(C4) 


The sum over j runs over all states and stands for the 
trace, I < j < 17, whereas the sum over i is also over 
all off-diagonal elements and corresponds to the matrix 
multiplication, thus 1 < i < 17^ = 289. The complicated 
indices of in the sum are a result of the vectorization 
of the trace. 

To estimate the ’leakage’ current ILeak mentioned in 
the main part of the paper we use the fact that the leak¬ 
age current does not depend on the state of the Ho atom 
and thus can be distinguished from Ixh- By compar¬ 
ing the height of the step of the differential conductance 
G = dI/dU{U) from elastic processes at zero bias 1/ = 0 
to voltages > AEsr at energies above the first excita¬ 
tion gap, the leakage current can be estimated. The step 
size is defined as s = G{U^)/G{Q) — 1. In the experiments 
a step of SExp ~ 0.9% was measured. Theoretically we 
expect without leakage, broadening and modulation as 
step size of srh ~ 10%. But if we assume that only 10% 
of the total current is described by Ixh, and 90% is leak¬ 
age current, we would expect a step size of s « 1%, which 
is in good agreement. For a more quantitative analysis of 
the dependence on the coupling constant ctb, and thus 
the current, we refer to Fig. |4]c). 
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